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We analyze further corrections due to the self-energy of the 
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(the pions); these corrections are of order a 5 In a for S states. 
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1 Introduction 



Exotic bound systems like pionium (the bound system of two oppositely charged pions) 
offer interesting possibilities for studies of fundamental properties of quantum mechanical bound 
states: the interplay between strong-interaction corrections and quantum electrodynamic cor- 
rections is of prime interest, and the small length scales characteristic of the heavy particles 
make it possible to explore effects of the virtual excitations of the quantum fields in previ- 
ously unexplored kinematical regimes |3|-|5). We do not wish to hide the fact that any potential 
high-precision experiments in this area are faced with various experimental difficulties. Our cal- 
culations address QED corrections to the spectrum of bound systems whose constituent particles 
are spinless; relativistic corrections to the decay lifetime of pionium have recently been discussed 
in H in the context of the DIRAC experiment at CERN. 

Here, we report on results regarding the spectrum of a bound system consisting of two spinless 
particles. We apply the simplified calculational scheme employed in J?J for the relativistic and 
recoil corrections to a bound systems of two "non-Dirac" particles to the case of two interacting 
spinless particles (see Sec. ^). We then recall known results on leading-order vacuum polariza- 
tion corrections in Sec. || and clarify the relative order-of-magnitude of the one- and two-loop 
electronic vacuum polarization, the relativistic and recoil corrections and the self-energy effects 
in pionium (also in Sec. |||). We then provide an estimate for the self-energy effect in Sec. ||, and 
we analyze the leading recoil correction of order a 5 (the Salpeter correction) which leads us to 
complete results for the scalar-QED logarithmic corrections of order a 5 In a. 



2 Breit Hamiltonian for Spinless Particles 

We start from the Lagrangian for a charged spinless field coupled to the electromagnetic field 
[see equations (6-50) - (6-51b) of g], 

C(x) = [(fy - ieAJ <f>*(x)} (0" + ieA») <j>{x) - m 2 f(x) 0(s) - J F^(x) F^(x) , (1) 

where the field strength tensor F^ u reads F fiv (x) = d fl A u (x) — dyA^x). We use natural Gaussian 
units with K = c = cq = 1. The transition current for a free spinless particle (A!* = 0) can be 
inferred from ((!]); it reads in momentum space 

f(p',p)=cj ) *(p')(pi»+pnHp)- (2) 

This current now has to be expressed in terms of nonrelativistic wave functions. Specifically, 
the j°-component has to reproduce the normalization of the nonrelativistic (Schrodinger) wave 
function. By contrast, according to Eq. (Q) the zero-component of the current reads 2m(j)*<j) in 
the nonrelativistic limit p'° — ► in ,p° — > in. The nonrelativistic wave functions are normalized 
according to 

' d d xcP* s (x)cP s (x) = l. (3) 



It is therefore evident that we cannot simply associate the relativistic wave function <f> with </>g; 
rather, we should define according to Eqs. (13) - (14) of 



V 2 p 

where p° = sjp 1 + m 2 ~ m is the energy of the free nonrelativistic particle (in deriving low- 
energy effective interactions, one always expands about /ree-particle amplitudes; all interactions 
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are treated as perturbations; note the analogy to nonrelativistic QED - NRQED - for spinor 
particles ||). The Klein-Gordon current, in the presence of external fields, reads in contrast to 

f(p',p) = <t>*{p') (p'^+p^ 1 -2eA») 0(p) . (5) 

The zero-component of this current can be interpreted as a charge density, which is not necessar- 
ily positive definite. Questions related to the normalization of the Klein-Gordon wave functions 
in this case are discussed in detail in Jl0|-13|. 

In terms of the Schrodinger wave function, the current is given as 



*s(pO&Cp): 



(6) 



f(p',p) 



p l + p" 
2 m 



4>s(p) , 



(7) 



where m is the mass of the particle. The atomic momenta p l and p H in Eq. (^) are of order 
Za. As shown below, interactions involving the spatial components f of the transition current 
give rise to relativistic contributions of order (Za) 4 to the spectrum. This is exactly the order 
of magnitude that is the subject of the current investigation. Therefore, although Eq. (0) is 
only valid up to corrections of relative order (Za) 2 , these can be neglected because the further 
corrections contribute to the energy levels at the order of (Za) 6 . Specifically, we can expect 
corrections proportional to (p l p 2 ) to the current f when a systematic expansion of the non- 
relativistic current is performed; these terms are analogous to those obtained for relativistic 
corrections to the current of spinor particles which can be obtained via a Foldy-Wouthuysen 
transformation [ 14 , . 

In the following, the index S on the wave function will be dropped, and the nonrelativistic am- 
plitudes describing the two interacting particles (with electric charges e% and e 2 ) will be denoted 
as <j)\ and </> 2 , respectively. Following the Breit Hamiltonian U(p 1: p 2 ,q) in momentum space 
is related to the invariant scattering amplitude M and to the photon propagator D^ v (q) in the 
following way [see also equation (83,8) in [fjl|]: 



M = e 1 e 2 j^(p' 1 ,p 1 )D^ u (q)j^(p' 2 ,p 2 ) 



ei e 2 



+ U(p 1 ,p 2 ,q) 



l(Pi)02(p 2 ) 



(8) 



where q = p' 2 — p 2 = —(Pi — P\)- We employ a Coulomb-gauge photon propagator, 



Doo(q) 



, D i:j (q) 



q 2 — ld 2 



5 ,J 



q l q> 



(9) 



where we can neglect the energy of the virtual photon for the derivation of next-to-leading order 
relativistic corrections, 



Dij(q) ~ 



S' J 



q l q> 



(10) 



The invariant scattering amplitude M then reads 



M 

d e2 



Wl)^(P2)^MPl)MP2) 



-^(pi)^(p' 2 ) 



Pi + Pi P2+P2 J_ 
2 mi 2 m 2 q 2 



5 ij 



1 l n3 



L(Pl)02(P2). (11) 



3 



We therefore identify 



U(p 1 ,p 2 ,q) 



e x e 2 (2p\ - q l ) (2p 3 2 + qi 



Anixm 2 



5 IJ 



e x e 2 f (2px - q) ■ (2p 2 + q) 



Am 1 7772 



(2 Pl • Q - q 2 ) (2p 2 -q + q 2 ) 



Px ■ P2 (Pi ' <?) (Pa • 9) 



eie 2 

77li 7772 

We now transform to the center-of-mass frame in which p 1 
for U(px,p 2 , q) becomes even simpler, 



P2 = P) so that the expression 



U(p, -p, q) 



ex e 2 

777 i 7772 



P 



ip-q) 



The formula (83,13) of [16| can now be employed in evaluating the Fourier transform, 



d 3 t 



(27T)* 



exp (lq ■ r) 



47r(a • q) (b ■ q) 



1 

2r 



a • b 



(a • r)(b • r) 



(13) 



(14) 



The Breit Hamiltonian, which we would like to denote by H^, is obtained by adding to the Fourier 
transform of ( |l3| ) the relativistic correction to the kinetic energy. Denoting with p = —id/dx 
the momentum operator in the coordinate-space representation, we obtain 



H B (r,p) 



8mf 



S 777o 



+ 



ex e2 



+ 



exe 2 r ■ [r ■ p) p 



8nr mx "^2 87rr 3 mxm 2 



(15) 



In the order of (Za) , there is no contribution due to virtual annihilation for spinless particles; 
corrections of this type would enter only for positronium and dimuonium || because they are 
caused by the spin-dependent part of the transition current [see Eqs. (83,20) and (82,22) of Ipf l, 
which is absent for spinless particles. For S states, virtual annihilation is altogether prohibited 
by angular momentum conservation. 

The matrix elements of the Breit Hamiltonian (15) for spinless particles can be evaluated on 
nonrelativistic bound states via computational techniques outlined in Sec. A3 of Ch. 1 of [[17]] . 
For mi = 1^2 = m and ex e 2 = —AirZa, we obtain 



Eb 



(Za) 2 m (Za) 4 m 



An 2 



2T7 3 



1 



21 + 1 



11 

3277 



(16) 



as the Breit energy for the energy levels of the bound system of two spinless particles, including 
relativistic corrections of order (Za) 4 . Here, we keep Z as a parameter which denotes the nuclear 
charge number in a bound system. Of course, for two particles each carrying an elementary 
charge, Z has to be set to unity. The fine-structure constant is denoted by a. The result ( |l6| ) 
agrees with previous calculations [f7, 18-^2|, notably with Eq. (38) of ipOj . 

It is instructive to compare the result ([!(]) with the known result for a single-particle system of 
mass m/2 satisfying the Klein-Gordon equation, bound to a nucleus with charge Ze. According 
to Eq. (2-86) of g], we obtain the "Klein-Gordon energy" (KG) 



Ekg 



(Za) 2 m (Za) 4 m 



An 2 



2 re 3 



1 

21 + 1 



.77 



(17) 
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The two results (16) and (17) are manifestly different in the order of (Za) 4 . 

^From (15) we conclude that the zitterbewegung term is absent for spinless particles. However, 



this statement is in need of further explanation because a considerable variety of physical inter- 
pretations exists in the literature with regard to the zitterbewegung term. We briefly expand: 
The Dirac a-matrices fulfills a = i[Hj),x] (Hd is the Dirac hamiltonian) as the relativistic 
generalization of the velocity operator. By contrast, in the nonrelativistic formalism, we have 
the analogous relation p/m = i[Hs,x] where H$ is the Schrodinger Hamiltonian. Since the en- 
matrices have eigenvalues ±1, the magnitude of the velocity of the electron - at face value - is 



equal to the velocity of light at any given instant. On p. 106 of |23 ], it is argued that "the expla- 
nation for this fact is that the electron carries out a fast irregular motion ( "zitterbewegung" ) - 
which is responsible for the spin - whereas the mean velocity is given by the momentum p/m". 
Note that the introduction of the Dirac matrix formalism is necessitated by the need to describe 
the internal degrees of freedom of the particle - the spin. On p. 71 of ||, it is shown that the 
zitterbewegung term can be traced to the positional fluctuations (5r 2 ) ~ 1/m 2 of the electron, 
and a connection is drawn to the Darwin term which results naturally in the context of the 
Foldy-Wouthuysen transformed Dirac hamiltonian. On pp. 117-118 of and p. 62 of ||, it is 
argued that the momentum p of a Dirac wave packet can be associated in a natural way with 
the group velocity, but that in addition to the group velocity term, there exist highly oscillatory 



terms which represent the zitterbewegung. Similarly, on pp. 139-140 of 24], it is shown that 
the zitterbewegung term can also be interpreted as arising from the interaction of the atomic 
electron with virtual electron-positron pairs created in the Coulomb field of the nucleus. This 
virtual electron-positron pair-creation is subject to the uncertainty principle and can occur only 
for time intervals of the order of At ~ 7i/(2mc 2 ) (where we temporarily restore the factor K). 
At the time the original atomic electron fills up the vacated negative-energy state (the bound- 
electron wave- function has negative-energy components), the escalated electron (which forms 
part of the virtual pair) is at most a distance cAt ~ 1/m away from the original electron. This 
distance is precisely of the order of magnitude of the fluctuations of the electron coordinate and 
consistent with the discussion on p. 71 of ||. All these interpretations elucidate different aspects 
of the same problem. 

In the context of the Breit hamiltonian, we would like to adhere to the definition that the zitter- 
bewegung term is the term of order (Za) 4 in the Breit Hamiltonian generated by a contribution 
which is manifestly proportional of 5{r) in coordinate space (or a constant in momentum space). 
Such a term is absent in the result (|i~5|). For spin-1/2 particles, such a term is generated by the 
multiplication of the photon propagator (proportional to 1/g 2 ) with the zero-component of the 
transition current which is given for a spin-1/2 particle as [see equation (4) of Q] 

ui'u = w 11-—^ + — — — \w. (18) 

\ 8m z im z J 

Here, u is the bispinor amplitude for the bound particle, and w is the bound-state Schrodinger 
wave function related by 

,2 



( 1--*. 

8 m' 



w 



(19) 

according to equation (3) of B. One might wonder why a term proportional to <5;o, apparently 



generated by a (5-function in coordinate space, prevails in the Breit energy (16). This term arises 
naturally when evaluating a matrix element of the structure (0g | (r ■ (r • p) p)/r 3 \rf>s} (last term 
of equation (|l5|) ) on the nonrelativistic wave function efts and should not be associated with the 
zitterbewegung. 
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Vacuum polarization corrections and self-energy effects, as well as corrections due to the strong 
interaction, are not included in (|l6l). These corrections will be discussed in the two following 
sections. 



3 Vacuum Polarization Effects 



As pointed out by various authors (e.g. [p^-|31[), the electronic vacuum polarization enters 
already at the order of a 3 [more precisely, a (Za) 2 ] in bound systems with spinless particles, 
because the spinless particles are much heavier than the electron, which means that the Bohr 
radius of the bound system is roughly of the same order of magnitude as the Compton wavelength 
of the electron. The Compton wavelength of the electron, however, is the fundamental length 
scale at which the charge of any bound particle is screened by the electronic vacuum polarization. 

The vacuum polarization (VP) correction to energy levels has been evaluated 



|g§|3-|g with 

nonrelativistic wave functions. We recall that the leading-order VP correction (due to the Uehling 
potential) can be expressed as 



AE = (ip\Vu\ip) = -C E E^, 

IT 



where 



E,; 



(Za) 



m 



An 2 



(20) 



(21) 



is the Schrodinger binding energy for a two-body system with two particles each of mass m 
[first term on the right-hand side of fll6l)1 • For the Ce coefficients, we recall the following known 
results @,[27|>f2§, 

Ce(1S) = 0.22, Ce(2S) = 0.10. (22) 

As an alternative to the nonrelativistic treatment, the Uehling correction could also be evaluated 
with relativistic bound-state Klein-Gordon wave functions (although these do not describe the 
two-body system accurately, as shown in Sec. ||). The difference could be interpreted as a rough 
estimate of further relativistic corrections not taken into account in the nonrelativistic treat- 
ment of the vacuum polarization. The result obtained by numerically solving the Klein-Gordon 
equation and numerically evaluating the Uehling correction is shown in the sixth row of Tab. |]. 
This relativistic result is in very good agreement with Eq. (p2|) . As pointed out in [32|, the strong 
interaction correction is also an a 3 effect (like the vacuum polarization) and enters at a relative 
order of a, i.e. on the level of about 1 % in pionium. 

We recall here also the known results for the vacuum polarization correction to the charge density 



at the origin |29| . |30|j 



A|fe(0)|- 
lV>is(0)| 2 



o 



1.36 and 



7T 



|V>2S(0)| 2 



a 



1.14. 



7T 



(23) 



Two-loop vacuum polarization effects enter at a relative order a 2 in pionium and are therefore 
of the same order of magnitude as the relativistic corrections mediated by the Breit interaction 
(discussed in Sec. ||| and termed "two-body correction" in Tab. |l|). The self-energy correction 
which is discussed in the following section is even smaller, but of considerable theoretical interest. 



4 Effects due to Scalar QED 



As shown in [33|, the leading logarithmic correction to the self energy can be obtained, in 



nonrelativistic approximation, from second-order perturbation theory based on nonrelativistic 
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quantum electrodynamics || (see also j34|). We will investigate here, in a systematic way, the 
leading logarithms generated for S states by self-energy and relativistic-recoil effects (the so- 
called Salpeter correction), and show that these are spin-independent. 

The quantized electromagnetic field is [see Eq. (5) of |3jJ], 



Mr) 



A=l,2 



d 3 k 



^(2ir) 3 2k 



e\(k) exp(— ifc • r) + A exp(ifc • r) 



(24) 



and the nonrelativistic interaction Hamiltonian for an atomic system with two spinless particles 
(charges e\ and e 2 and masses m\ and 7712) reads 



Hi = -—pi ■ A(n) + ^A(ri) 2 
mi Irrii 



— P2 ■ A(r 2 ) + -^-A(r 2 ) 5 
m 2 2m 2 



For two spin-1/2 particles, the terms 

_£]_ 

mi 



<Ti ■ B(ri) 



e 2 
m 2 



a 2 • B(r 2 ) 



(25) 



(26) 



have to be added to Hi [see Eq. (7) of |p3|| ]. We will carry out the calculations for the general 
case of one particle of charge e\ = e and the other having a charge e 2 = — Ze (we follow the 
convention of || that for hydrogen, e is the physical charge of the electron, i.e. e = — |e|). The 
unperturbed Hamiltonian of the system of the two particles and the electromagnetic field reads 



[see e.g. Eq. (6) of ||], 
Ho 



pi 



+ 



2mi 2m 2 



Za \ ^ 

A=l,2 



drkka 



k,x a k,\ 



(27) 



where r = r\ — r 2 . The eigenstates of the "atomic part" H$ of this Hamiltonian in the center- 
of-mass system pi + p 2 = are the nonrelativistic Schrodinger-Coulomb wave functions for a 
reduced mass m r = mim 2 /(rrti + m 2 ) [here, the "atomic part" Hq excludes the photon field, 
i.e. the last term of (\ZT\j\. We denote p = p 1 = —p 2 - 

Given that the first-order perturbation {(f)s\Hi\(f)s) vanishes, the second-order perturbation yields 
the dominant nonvanshing perturbation. When evaluated on an atomic state, it is given by 



SE SE 



k\Hi 



Ho-E s 



Hilfc). 



(28) 



The interaction Hamiltonian ( |25| ) gives rise to QED corrections that involve both particles (in 
the current context, these are recoil corrections involving the product e\ e 2 ), and also to terms 
which involve only a single particle and are proportional to e\ or e 2 . The latter effects correspond 
to the self-energies of the two particles. 

The low-energy part of the self-energy in leading order p5l can be inferred directly from fl2§|), 
and it can be seen that the spin-dependent parts from ( P6[) vanish in the leading order in the 
(Za)-expansion [33|: 



6tt 2 



dkk 



P 



P 



mi Hq 4 - - (E s - k) mi 



+ (ei <-> e 2 ,mi 



m 2 



where 



E s 



(Za) 2 m T 
2^ 



(29) 



(30) 
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is the Schrodinger energy (m T in the reduced mass of the atomic system under investigation). 
Starting from the spin-independent expression (^9|), it is now relatively straightforward to show 
that the leading "self-energy logarithm" for S states is given by 



-2 



$E S E ~ — 5 3 — SlQ 

6 Tin 6 



,3 ™3 



a{Zaf-^ + Z{Zaf 

mi ms 



(31) 



This result is by consequence spin-independent. The derivation is simplified when using the 
e-method developed and used in various bound-state calculations |14,35|, 36j. The two terms 



in square brackets in (plD correspond to the two self-energies of the two constituent particles 
with charges e\ = e and ei = —Ze and masses m\ and rri2, respectively. It has been pointed 
out |?7) that in contrast to the self-energy corrections, the vacuum polarization corrections given 
in Eq. ( |23"| ) must not be double-counted. The "double-counting" of self-energy corrections (and 
lack of it in the vacuum-polarization case) finds a natural explanation in our formalism: whereas 
the vacuum-polarization correction mainly leads to a modification of the 1/r-type Coulomb at- 
traction in ( |27|) within a nonrelativistic effective theory, the structure of the interaction Hamil- 
tonian ( |25| ) implies the existence of the two self-energies of the two constituent particles of the 
atomic system. 

It might be instructive to point out that the formula ([H]) is consistent with Welton's argument 
for estimating the self-energy effect on a bound particle which is based on analyzing the influence 
of the fluctuating electromagnetic field [a detailed discussion is given on pp. 80-82 of ||]. For a 
system with two particles of equal mass m\ = mi = m, we have m T = m/2. 

The leading-order recoil correction (Salpeter correction) can also be inferred from the interaction 
Hamiltonian (^) via second-order perturbation theory, by "picking up" terms that involve 
products e\ ei. It has been shown in |33| that the leading logarithm (for S states) of the Salpeter 
correction is spin- independent (just like the leading logarithm of the self-energy correction). The 
Salpeter correction is usually referred to as a relativistic recoil (RR) correction. By following [^] , 
we obtain for the leading logarithm of this effect 

2{Zaf ( 1 \ m 3 r 

$E RR » 3- 610 In — . (32) 

Snn \Za J m\mi 

This correction involves only the products e\ ei = —AirZa and can therefore be written as a 
function of Za alone. 

For pionium, we have Z = 1, m\ = mi = m = m n , m T = m^/2. The leading logarithmic correc- 
tion from scalar QED for pionium in the order of a 5 In a is obtained by adding the corrections 
(|3lD and (||), 

3 a 5 / 1 \ 

5E\ og = 5E S e + SE RR = 5- In - m T . (33) 

4 7r n 6 \Oi J 

The non-logarithmic term of order a 5 is spin-dependent, and its evaluation requires a relativistic 
treatment of the self-energy effect of a bound spinless particle; such a calculation would be of 
considerable theoretical interest, but the size of the effect for pionium, which is roughly two 
orders of a smaller than the leading vacuum polarization correction, precludes experimental 
verification in the near future. However, we would like to point out here that a fully relativistic 
treatment of this problem, including a detailed discussion of the renormalization of the self 
energy of the spinless particle, has not yet been accomplished. Scalar QED is a renormalizable 
theory g. 

The dominance of vacuum polarization over self-energy effects in pionium is expressed, in partic- 
ular, by the fact that even two-loop vacuum polarization of order a 4 has a stronger effect on the 
spectrum of pionium than the leading logarithm from Eq. (^), and that the strong-interaction 
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correction of order a 3 [32] has to be well understood before any experimental verification of ( |33| ) 
appears feasible. Finally, we remark that for a manifestly non-elementary particle like the pion 
which has a finite charge radius, form-factor corrections have to be taken into account. 



Table 1: QED contributions to the IS level of pionium in eV. For a 
further discussion of the corrections see the text. 



Contribution 


Energy (eV) 


One-body Klein-Gordon [Eq. (17)] 


-1858.19895 


Higher Order Klein-Gordon [Exact - Eq. fll7|)] 


-0.00001 


Form factor correction to Klein-Gordon (0.61 fm) 


0.01308 


Two-body correction (Breit) [Eq. ([i~6|) - Eq. (|i~7|)1 


0.04329 


Uehling (with a relativistic wave function) 


-0.94228 


Vac. Pol. (Wichmann-Kroll) 


0.00001 


Vac. Pol. (Kallen-Sabry) 


-0.00729 


Vac. Pol. (iterated Uehling) 


-0.00113 


Self-Energy [Eq. @)] 


0.00302 


Salpeter correction [Eq. (|33[)] 


0.00038 


Total 


-1859.08986 



5 Numerical evaluation of QED corrections 

In order to provide a more complete picture of pionium we have evaluated numerically a number 
of QED corrections to the IS level of pionium. We explicitly exclude QCD corrections whose 



evaluation represents a difficult separate problem p2 |. We proceed as follows (see the sequence 
of the rows of Tab. |l|): 

• We start from the one-body Klein-Gordon energy for a particle of (reduced) mass m/2, 
given in Eq. (|j"7|), including relativistic corrections of order (Za) 4 . 

• We use a Klein-Gordon equation numerical solver developed for pionic atoms which was 
developed originally for the precise evaluation of vacuum polarization corrections [38|, 36] 



in order to supplement the (almost neglible) difference between the exact one-body rela- 
tivistic Klein-Gordon energy and the (Za) 4 -result from Eq. ([17]), thereby confirming the 
expression ( |i~7l) for the relativistic correction. 



• The effect due to the pion Coulomb form factor is also included in an approximative frame- 
work by replacing in the numerical solution of the Klein-Gordon equation the Coulomb 
potential by the interaction potential of two uniformly charged spheres of mean-square 
radius R rms = 0.61 fm. For the calculation, we employ the radius and the pion mass from 
the particle data grou p |40| ] (all other physical constants used in the evaluations come from 
the 1998 adjustment gf^ 

• We add as a "two-body correction" the difference of the results from Eq. di~6| ) for the 
energy of the relativistic two-body system and the one-body result given in Eq. (|i~7"l). 
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The Uehling potential is evaluated in a relativistic framework, using the relativistic numer- 
ical equation solvers |38,3S]. The result is in very good agreement with the nonrelativistic 



treatment discussed in Sec. || (this is natural because Z = 1) 

The higher-order VP corrections attributed to Wichmann-Kroll p2| and Kallen-Sabry 
are supplemented, as well as an evaluation of the iterated loop-after-loop Uehling contri- 
bution to all orders in a. The Wichmann-Kroll correction is here negligible because 



Z = 1. VP potentials given in Ref. [44] and analytic expressions from [45| are used 



More details about the numerical procedure to evaluate these corrections can be found in 



Ref. H| 



Finally, the effects due to the scalar self-energy given in Eq. (ph) and the relativistic recoil 



(Salpeter) correction listed in Eq. (|33| ) are added. 



The main observation one can draw from Tab. [l| is that the vacuum-polarization is 300 times 
larger than the self-energy. Even the Kallen and Sabry correction, while of order a? is about 
twice as large as the self-energy correction. The iterated Uehling correction, while dominated by 
terms of order a 2 is dominated in turn by the scalar self-energy. Nevertheless, we stress here that 
the numbers contained in Tab. |l] will be modified when the Klein-Gordon equation is solved with 
the strong interaction potential incorporated directly into the equation solver. Analogously, one 
cannot avoid having to solve the Dirac equation exactly for high-Z systems where the electron 
wave function rests significantly inside the nucleus piffl . At present, in view of the formidable 
experimental difficulties associated with a study of the atomic spectrum of pionium, we give the 
numbers in Tab. [I] as an indication of the relative size and order-of-magnitude of the specific 
QED corrections. 



6 Conclusion 

We have presented in Sec. ^ a simplified derivation for the relativistic and recoil corrections of or- 
der a 4 to a bound state of two spinless particles. The results agree with previous calculations [p(|. 
As evident from equation (|l5l) and discussed in further detail in Sec. the zitterbewegung term 
is absent in a bound system of two spinless particles. 

The self-energy effect is suppressed in systems with spinless particles in comparison to the 
vacuum polarization effect as discussed in Sees. |3| and |], because the lightest known spinless 
particle is much heavier than the electron, which implies that the electronic vacuum polarization 
effect is larger by two orders of a than the self-energy effect in bound systems of spinless particles. 
We provide a complete result for the leading scalar-QED correction in pionium of order a 5 In a 
in Sec. ^. A list of further QED corrections to the IS level of pionium is presented in Tab. |l|. 
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